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Abstract We investigate the problem of designing a minimum 
cost flow network interconnecting n sources and a single sink, each with 
known locations and flows. The network may contain other unprescribed 
nodes, known as Steiner points. For concave increasing cost functions, 
a minimum cost network of this sort has a tree topology, and hence can 
be called a Minimum Gilbert Arborescence (MGA). We characterise the 
local topological structure of Steiner points in MGAs, showing, in par- 
ticular, that for a wide range of metrics and cost-functions the degree of 
each Steiner point is 3. 

Keywords: Gilbert network, minimum cost networks, network flows, 
Steiner trees 

1 Introduction 

The Steiner Minimum Tree (SMT) problem asks for a shortest network 
spanning a given set of nodes {terminals) in a given metric space. It differs 
from the minimum spanning tree problem in that additional nodes, referred 
to as Steiner points, can be included to create a spanning network that is 
shorter than would otherwise be possible. This is a fundamental problem 
in physical network design optimisation, and has numerous applications, 
including the design of telecommunications or transport networks for the 
problem in the Euclidean plane (the h metric), and the design of microchips 
for the problem in the rectilinear plane (the l\ metric) j5]. 

Gilbert [I] proposed a generalisation of the SMT problem whereby sym- 
metric non-negative flows are assigned between each pair of terminals. The 
aim is to find a least cost network interconnecting the terminals, where each 
edge has an associated total flow such that the flow conditions between ter- 
minals are satisfied, and Steiner points satisfy Kirchhoff's rule (ie, the net 
incoming and outgoing flows at each Steiner point are equal). The cost of 
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an edge is its length multiplied by a non-negative weight. The weight is 
determined by a given function of the total flow being routed through that 
edge, where the function satisfies a number of conditions. The Gilbert net- 
work problem (GNP) asks for a minimum-cost network spanning a given set 
of terminals with given flow demands and a given weight function. 

A variation on this problem that we will show to be a special case of 
the GNP occurs when the terminals consist of n sources and a unique sink, 
and all flows not between a source and the sink are zero. This problem is of 
intrinsic interest as a natural restriction of the GNP; it is also of interest for 
its many applications to areas such as drainage networks [7] , gas pipelines p] , 
and underground mining networks [2]. 

If the weight function is concave and increasing, the resulting minimum 
network has a tree topology, and provides a directed path from each source 
to the sink. Such a network can be called an arborescence, and we refer to 
this special case of the GNP as the Gilbert arborescence problem (GAP). 
Traditionally, the term 'arborescence' has been used to describe a rooted 
tree providing directed paths from the unique root (source) to a given set 
of sinks. Here we are interested in the case where the flow directions are 
reversed, i.e. flow is from n sources to a unique sink. It is clear, however, 
that the resulting weights for the two problems are equivalent, hence we will 
continue to use the term 'arborescence' for the latter case. Moreover, if we 
take the sum of these two cases, and rescale the flows (dividing flows in each 
direction by 2), then again the weights for the total flow on each edge are 
the same as in the previous two cases. This justifies our claim that the GAP 
can be treated as a special case of the GNP. It will be convenient, however, 
for the remainder of this paper to think of arborescences as networks with 
a unique sink. 

A minimum Gilbert arborescence (MGA) is a (global) minimum-cost ar- 
borescence for a given set of terminals and flows, and a given cost function. 
All flows in the network are directed towards the unique sink. In this paper 
we investigate the local topological structure of Steiner points in MGAs, 
over a wide range of different metrics and cost-functions. Such an under- 
standing of the topological structure is a key step towards designing efficient 
algorithms for constructing MGAs. 

In Section 2 we specify the nature of the weight function that we consider 
in this paper, and formally define minimum Gilbert networks and Gilbert 
arborescences in Minkowski spaces (which generalise Euclidean spaces). In 
Section 3 we give a general topological characterisation of Steiner points in 
such networks, for smooth Minkowski spaces. We then apply this charac- 
terisation, in Section 4, to the smooth Minkowski plane with a linear weight 
function to show that in this case all Steiner points have degree 3. In Sec- 
tion 5 we derive a similar result in higher dimensional Euclidean spaces for 
a wide class of weight functions. 
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2 Preliminaries 



2.1 Minkowski spaces and Steiner trees 

The cost functions for the networks we consider in this paper make use of 
more general norms than simply the Euclidean norm. Hence, we introduce a 
generalisation of Euclidean spaces, namely finite-dimensional normed spaces 
or Minkowski spaces. See j!2| for an introduction to Minkowski geometry. 

A Minkowski space (or finite- dimensional Banach space) is M n endowed 
with a norm ||-||, which is a function ||-|| : M n — > R that satisfies 

• H^ll > for all x € M n , ||x|| = only if x = 0, 

• \\ax\\ = |a|||x|| for all a £ M and x £ M. n , and 

• \\x + y\\ < \\x\\ + ||y||. 

We use ||-||2 to denote the Euclidean (I2) norm. 

We now discuss some aspects of the SMT problem, since this is a special 
case of the GNP, where all flow are zero. Our terminology for the SMT 
problem is based on that used in [5|. Let T be a network interconnecting 
a set N = {pi, . . . ,p n } of points, called terminals, in a Minkowski space. 
A vertex in T which is not a terminal is called a Steiner point. Let G{T) 
denote the topology of T, i.e. G{T) represents the graph structure of T but 
not the embedding of the Steiner points. Then G{T) for a shortest network 
T is necessarily a tree, since if a cycle exists, the length of T can be reduced 
by deleting an edge in the cycle. A network with a tree topology is called 
a tree, its links are called edges, and its nodes are called vertices. An edge 
connecting two vertices a, b in T is denoted by ab, and its length by \\a — b\\. 

The splitting of a vertex is the operation of disconnecting two edges 
av, bv from a vertex v and connecting a, b, v to a newly created Steiner point. 
Furthermore, though the positions of terminals are fixed, Steiner points can 
be subjected to arbitrarily small movements provided the resulting network 
is still connected. Such movements are called perturbations, and are useful 
for examining whether the length of a network is minimal. 

A Steiner tree (ST) is a tree whose length cannot be shortened by a 
small perturbation of its Steiner points, even when splitting is allowed. By 
convexity, an ST is a minimum-length tree for its given topology. A Steiner 
minimum tree (SMT) is a shortest tree among all STs. For many Minkowski 
spaces bounds are known for the maximum possible degree of a Steiner point 
in an ST, giving useful restrictions on the possible topology of an SMT. For 
example, in Euclidean space of any dimension every Steiner point in an ST 
has degree three. Given a set N of terminals, the Steiner problem (or Steiner 
Minimum Tree problem) asks for an SMT spanning N. 

2.2 Gilbert flows 

Gilbert [1] proposed the following generalisation of the Steiner problem in 
Euclidean space, which we now extend to Minkowski space. Let T be a net- 
work interconnecting a set N = {p\, . . . ,p n } of n terminals in a Minkowski 
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Figure 1: An example where split-routing is cheaper 

space. For each pair Pi,Pj, i ^ j of terminals, a non-negative flow t- L j = tji 
is given. The cost of an edge e in T is w(t e )l e , where l e is the length of e, 
t e is the total flow being routed through e, and w(-) is a unit cost weight 
function defined on [0, oo) satisfying 

w(0) > and w(t) > for all t > 0, (1) 

w(t 2 ) > w(h) for all t 2 > h > 0, (2) 

w{h + t 2 ) < w(ti) +w{t 2 ) for all h,t 2 >0, (3) 

w(-) is a concave function. (4) 

That the function w is concave means by definition that — w is convex. A 
network satisfying Conditions Q^(|3| is called a Gilbert network. For a given 
edge e in T, w(i e ) is called the weight of e, and is also denoted simply by 
w e . The fotaZ cost of a Gilbert network T is the sum of all edge costs, i.e. 

C(T) = £>(i e )Ze 

where £ is the set of all edges in T. A Gilbert network T is a minimum 
Gilbert network (MGN), if T has the minimum cost of all Gilbert networks 
spanning the same point set N, with the same flow demands tij and the 
same cost function w(-). By the arguments of [3], an MGN always exists 
in a Minkowski space when Conditions are assumed for the weight 

function. 

Conditions ([TJ-Q above ensure that the weight function is non-negative, 
non-decreasing and triangular, respectively. These are natural conditions 
for most applications. Unfortunately, the first three conditions alone do not 
guarantee that a minimum Gilbert network is a tree. To show this, we now 
give an example of a Gilbert network problem with two sources and one sink 
in the Euclidean plane, where there exists a split-route flow that has a lower 
cost than any arborescence. 

For this example there are two sources pi,p 2 and a sink q which are 
the vertices of a triangle Apip 2 q with edge lengths \\p\ — p 2 \\ 2 = 1 and 
||Pl — 1 1 2 = \\P2 — ?|b = 10, as illustrated in Figure [TJ The tonnages at p\ 
and p 2 are 2 and 4, respectively. The weight function is w(t) = \(3t + l)/2] , 
i.e., (3i + l)/2 rounded up to the nearest integer. This function is positive, 
non-decreasing and triangular, but not concave. For the example we only 
need the following values: 
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Figure 2: The minimum Gilbert arborescence 




Figure 3: Constructing the weighted Fermat-Torricelli point 

Routing 1 unit of the flow from p2 via p\ to q gives a Gilbert network 
(Figure [T]) of total cost 

w(l)\\pi -P2II2 + w(3)\\pi - q\\ 2 + w(3)\\p 2 - qh = 102 - 

For a Gilbert arborescence we route the tonnages from p\ and pi to q via 
some point s as in Figure I) We calculate a minimum Gilbert arborescence 
by using the weighted Melzak algorithm as described in [1]. To construct 
the weighted Fermat-Torricelli point s, first construct the unique point p 
outside Apip2q such that \\p — P1W2 = 0.7 and \\p — P2II2 = 0-4 (Figure [3]). 
Then construct the circumscribed circle of App±p2, which will intersect the 
so-called weighted Simpson line pq in the required point s. The resulting 
total cost is 

w(2)\\pi - s|| 2 + io(4)||p2 - s\\2 + w(6)\\s - q\\ 2 

= w(6)\\p-q\\ 2 = \/9982.5 + 7\/3890.25 
= 102.074 

as some trigonometry will show. This shows that split routing can be nec- 
essary when the weight function is not concave. 

For the remainder of the paper we assume that the weight function w is 
concave in addition to the other three conditions. In this case it is known 
jU [H] that in the case where there is a single sink there always exists a 



5 



minimum Gilbert network that is a Gilbert arborescence. This means that 
we can (and will) without loss of generality only consider MGAs. (Note that 
in j3], Condition ([3]), which we call the triangular condition, was incorrectly 
interpreted as concavity of the cost function.) 

The Gilbert network problem (GNP) is to find an MGN for a given ter- 
minal set N, flows tij and cost function w(-). Since its introduction in 
various aspects of the GNP have been studied, although the emphasis has 
been on discovering geometric properties of MGNs (see [3], [11], [13], pH]). 
As in the Steiner problem, additional vertices can be added to create a 
Gilbert network whose cost is less than would otherwise be possible, and 
these additional points are again called Steiner points. A Steiner point s in 
T is called locally minimal if a perturbation of s does not reduce the cost 
of T. A Gilbert network is called locally minimal if no perturbation of the 
Steiner points reduces the cost of T. 

The special case of the Gilbert model that is of interest in this work 
is when N = {pi, . . . ,p n , q} is a set of terminals in a Minkowski space, 
where pi, . . . ,p n are sources with respective positive flows t\, . . . ,t n , and 
q is the sink. All flows are between the sources and the sink; there are 
no flows between sources. It has been shown in [llj that concavity of the 
weight function implies that an MGN of this sort is a tree. Hence we refer 
to an MGN with this flow structure as a minimum Gilbert arborescence 
(MGA), and, as mentioned in the introduction, we refer to the problem of 
constructing such an MGA as the Gilbert arborescence problem (GAP). 

If vi and V2 are two adjacent vertices in a Gilbert arborescence, and flow 
is from v\ to Vi then we denote the edge connecting the two vertices by v\v<i- 

3 Characterisation of Steiner Points 

In this section, we generalise a theorem of Lawlor and Morgan |6 to give a 
local characterisation of Steiner points in an MGA. The characterisation in 
[6] holds for SMTs, which correspond to the case of MGAs with a constant 
weight function. Their theorem is formulated for arbitrary Minkowski spaces 
with differentiable norm. Our proof is based on the proof of Lawlor and 
Morgan's theorem given in [8]. A generalisation to non-smooth norms is 
contained in [TO] for SMTs and in [15] for MGAs. Such a generalisation is 
much more complicated and involves the use of the subdifferential calculus. 

We first introduce some necessary definitions relating to Minkowski ge- 
ometry, in particular with relation to dual spaces. For more details, see 



We denote the inner product of two vectors x, y E W 1 by (x, y). For any 
given norm ||-||, the dual norm ||-||* is defined as follows: 




We say that a Minkowski space {W 1 , ||-||) is smooth if the norm is differ- 
entiable at any x ^ o, i.e., if 




sup (z, x) . 



lim 



x + th\\ — \\x 



: Uh) 



t 
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Figure 4: A Gilbert network with star topology. 



exists for all x, h € M n with x ^ o. It follows easily that f x is a linear 
operator f x : W 1 — > M and so can be represented by a vector x* G M" - , called 
the dual vector of x, such that (x*,y) = f x (y) for all y € W 1 , and ||x* ||* = 1. 
In fact x* is just the gradient of the norm at x, i.e., x* = V||x||. 

More generally, even if the norm is not differentiable at x, a vector 
x* £ W 1 is a dual vector of x if x* satisfies (x*,x) = \\x\\ and ||x*||* = 1. By 
the Hahn-Banach separation theorem, each non-zero vector in a Minkowski 
space has at least one dual vector. A Minkowski space is then smooth if and 
only if each non-zero vector has a unique dual vector. 

A norm is strictly convex if ||x|| = \\y\\ = 1 and x ^ y imply that 
\\n(x + y)\\ < 1, or equivalently, that the unit sphere 

5(||-||) = {i£l n : ||x|| =1} 

does not contain any straight line segment. A norm ||-|| is smooth [strictly 
convex] if and only if the dual norm ||-||* is strictly convex [smooth, respec- 
tively] . 

Theorem 1. Suppose a smooth Minkowski space (W 1 , ||-||) is given together 
with a weight function w that satsfies Conditions ([TJ-Q, sources pi, ... ,p n G 
W 1 , and a single sink q £ M™, all different from the origin o. Let the flow 
associated with pi be tj. (See Figure Q ) For each p\ let p* denote its dual 
vector, and let q* denote the dual vector of q. Then the Gilbert arborescence 
with edges opi, i = 1, . . . , n and oq, where all flows are routed via the Steiner 
point o, is a minimal Gilbert arborescence if and only if 

n n 

^w(U)Pi+w(jr,ti)q* = o (5) 

i=l i=l 

and 

WY^wiUMW* < w(J2ti) for all I <Z {1, . . . ,n}. (6) 

iei iei 

Note: We think of Condition [5] as a flow-balancing condition at the 
Steiner point, and Condition [6] as a condition that ensures that the Steiner 
point does not split. 
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Proof. (=>) We are given that the star is not more expensive than any other 
Gilbert network with the same sources, sink, flows and weight function. 

In particular, o is the so-called weighted Fermat-Torricelli point of the 
n + 1 points pi, . . . ,p n , q with weights t\,...,t n , Y17=i respectively, which 
implies the balancing condition ([5]). We include a self-contained proof for 
completeness. If the Steiner point o is moved to — te, where iGl and e E M. n 
is a unit vector (in the norm), the resulting arborescence is not better, by 
the assumption of minimality. Therefore, the function 



¥e(t) = ^ W ( t i)(.\\Pi + te \\ ~ \\Pi 
i=l 

n 



(5>)(ll9 + *e||-|M|)>0 



w , 

=1 



attains its minimum at t = 0. For t in a sufficiently small neighbourhood of 
0, pi + te ^ o and q + te / o, hence <p e is differentiable. Therefore, 



o = ip'M 




Pi + te\\ - \\p,j 
t 



i=i 



\q + te\\ 
t 



J2w(ti){p*,e) + w(J2u) (q*,e) 

8=1 1=1 
In n \ 



,1=1 1=1 



Since this holds for all unit vectors e, §5§ follows. 

To show (|6| for each / C {1, . . . ,n}, we may assume without loss of 
generality that 1^0 and / ^ { 1 , . . . , n} . Consider the Gilbert network 
obtained by splitting the Steiner point into two points o and +te (t £ 1, e a 
unit vector) as follows. Each pi, i £ I, is still adjacent to o with flow ij, and 
q is joined to o with flow Y17=i but now each pi, i £ I, is adjacent to te 
with flow tj, and te is adjacent to o with flow ^ieJ" ^> as shown in Figure [5j 
Since the new network cannot be better than the original star, we obtain 
that for any unit vector e, the function 

n 

^e(i) = ^2w(ti)(\\pi - te\\ - ||pi||) + to(5^*i)|t| > 

i£l i=l 

attained its minimum at t = 0. Although ip e is not differentiable at 0, we 
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Figure 5: The Gilbert network obtained by splitting the Steiner point o. 



can still calculate as follows: 
< lim 



\\pi - te\\ - \\pj\ 
t 

iei i=i 



t 5g|j.E li, (**)~ ~ t —+ w (%2 t i) 



I i=i 

Therefore, (Yliel w (^i)Pi j e ) — ^(EiLi *») ^ or an un ^ vectors e, and ([6| 
follows from the definition of the dual norm. 

(<=) Now assume that p\,p* n ,q are dual unit vectors that satisfy ^ and 
Q. Consider an arbitrary Gilbert arborescence T for the given data. For 
each i, let Pj be the path in T from p, to q, i.e., Pj = ij 4' ■ • • where 

xf' = P i, x^ = q, and XjXj^ are distinct edges of T for j = 1, . . . , ki — 1. 
For each edge e of T, let S e = {i : e is on path Pj}. Then the flow on e is 
X^eS e and the total cost of T is 

w (^2 t i)\\x-y\\. 

e=xy is iG5 e 
an edge of T 

The cost of the star is 

n n 

Y,w{ti)\\Pi\\+w{Y,ti)\\q\\ 
i=i i=i 

n n 

= J2w(ti) (p*,pi) + wC^U) (q*,q) 



8=1 1=1 



8=1 

n fcj— 1 

(ft*'*? 

8=1 J=l 
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e=xy is \ieS e I 
an edge of T 

e=xy is «s5 e 
an edge of T 

e=xy is iGS e 
an edge of T 

Note that the necessity of the conditions ^ and ([6| holds even if the 
weight function is not concave. It is only in the proof of the sufficiency that 
we need all minimal Gilbert networks with a single source to be arbores- 
cences. 



4 Degree of Steiner Points in a Minkowski plane 
with linear weight function 

We now apply the characterisation of the previous section in the two-dimen- 
sional case, assuming further that the weight function is linear: w(t) = d+ht, 
d>0,h>0. 

Theorem 2. In a smooth Minkowski plane and assuming a linear weight 
function w(t) = d + ht, d > 0, h > 0, a Steiner point in an MGA necessarily 
has degree 3. 

Proof. By Theorem [TJ an MGA with a Steiner point of degree n + 1 exists 
in IR 2 with a smooth norm ||-|| if and only if there exist dual unit vectors 
p*, . . . q* G R 2 such that 

n n 

^(d + hti)p* + (d + hY^ U)q* = o 
i=i i=i 

\\^2{d + hti)p*\\* Kd + hJ^ti for all/C {l,...,n}. 

iei iei 



and 



Label the p* so that they are in order around the dual unit circle. Let 
v* = (d + hti)p* and w* = (d + h Y12=l U)q* ■ Then the conditions become 

v*-\ h v* n + w* = o, 

and 

< d + hj^ti for all JC{1 n}. (7) 

iei iei 

Thus we may think of the vectors , . . . , , w* as the edges of a convex 
polygon with vertices a*j = Yjl=i v l j 3 = 0> • • • > n m this order (see Figure 6). 

Assume for the purpose of finding a contradiction that n > 3. Then 
the polygon has at least 4 sides. Note that the diagonals a^a* and a*_ 1 a* 
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Figure 6: A polygon with edges corresponding to v*. 

intersect. Applying the triangle inequality to the two triangles formed by 
these diagonals and the two edges v* and w* (as illustrated in Figure |6j), we 
obtain 

\\dj\\* + \\o* n — OjLiH* > \\v*j\\* + ||w*||* 

n 

= d + htj + d + hy^Jj 

i=l 

j n 

= d + hy^U + d + 

i=l i=j 

>nE<ir + iE<ir b y @ 

i=l i=j 
= HIP + IK - a *j-l\\*- 

Therefore, equality holds throughout, and we obtain equality in the triangle 
inequality. Since the dual norm is strictly convex, it follows that v*- and 
w* are parallel. This holds for all j = 2, . . . , n — 1. It follows that p\ = 
■ ■ ■ = Pn = ~Q* ■ Geometrically this means that the unit vectors ]]^T||Pi an d 
— jj^jj-g all have the same supporting line on the unit ball. We can think of 
this condition on the vectors pi and q as a generalisation of collinearity to 
Minkowski space. 

Choose a point S2 on the edge opi such that the line through S2 parallel 
to op n intersects the edge op\ in si, say, with s\ ^ o. See Figure [7j Because 
of the straight line segments on the boundary of the unit ball, \\x + y\\ = 
\\x\\ + \\y\\ for any x, y such that the unit vectors jrar^ and w^ny lie on this 
segment. In particular, 

H-S2 — + [| Si — o|| = ||S2 — 0||. (8) 

Now replace P20 by the edges P2S2 and S2S1 , replace p\o by p\S\ and sio, and 
add the flow £2 to s\o. The change in cost in the new Gilbert arborescence 
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Figure 7: Illustration of the proof of Theorem [2] for a unit ball with straight 
line segments on the boundary. 



is 

(™(*i)lbi - sill + w{h + t 2 )\\si - o\\ + w{t 2 )\\p 2 - s 2 \\ + w(t 2 )\\s 2 - Si||) 
- (w(h)\\pi - o|| - w(t 2 )\\p 2 - o\\) 
= -«;(ti)[|ai[| -w(i 2 )||s 2 || +w(t! +t 2 )\\ Sl \\ + w(t 2 )(\\s 2 \\ - \\ Sl \\) by @ 

= {w{h + 1 2 ) - w{h) - w{t 2 ))\\ Sl \\ 

= {d + h(h + t 2 ) - (d + hti) - (d + ht 2 ))\\si\\ 
= -d\\sx\\ < 0. 

We have shown that a Steiner point of degree at least 4 leads to a decrease 
in cost. In an MGA a Steiner point must then necessarily be of degree 3. □ 



5 Degree of Steiner points in Euclidean space 

We now consider Gilbert arborescences in Euclidean space (of arbitrary di- 
mension) with the weight function w(t) = d + ht a , where d, h > and 
< a < 1. Note that for all these values of a, the weight function w satis- 
fies all of the conditions (JTJ-Q. We show that if < a < 1/2 or a = 1, then 
the maximum degree of a Steiner point is 3, while for each a S (1/2, 1) we 
provide an example of an MGA with a Steiner point of degree 4. This exam- 
ple is three-dimensional, and the tonnage goes to infinity as a approaches 
1/2 or 1. When a = 1 then the previous section shows that Steiner points 
are necessarily of degree 3 in the Euclidean plane. Although we have no 
examples of degree 4 Steiner points in the Euclidean plane for other values 
of a, we consider their existence to be highly unlikely. 

We show that there is a very general class of weight functions for which 
the Steiner points are necessarily of degree 3 (Theorem [4] below). Our proof 
is completely independent of the dimension. 

As previously, we denote the Euclidean norm by ||-||2, the inner product 
of two vectors u and v by (u, v), and the weight function by w. 

We first reformulate Theorem [T] for Euclidean space. 

Corollary 3. A Steiner point of degree m + 1 is possible in some MGA in 
Euclidean space with a weight function w(-) that satisfies Conditions ([T|- 
Q, if and only if there exist vectors v\, . . . , v m and tonnages t\, . . . , t m > 
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such that 



IKH2 = w(U), i = l,...,m (9) 

m m 

lE«i||2 = (10) 

i=l i=l 

V7C {l,...,m} with 2 < \I\ < m - 2, Vi|| 2 < t»). (11) 

Theorem 4. If the weight function w( ) satisfies Conditions (JIJ-Q, and 
is differentiable with (w 2 )' convex and w(0) > 0, then all Steiner points in 
MGAs have degree 3. 

Note that the hypothesis is indeed satisfied for the weight function w(t) = 
d + ht a for any d, h > and a £ [0, 1/2] U {1}, but not when a £ (1/2, 1). 

For the proof of Theorem [4] we need the following inequality valid for 
functions with convex derivative. 

Lemma 5. If f : [0, 00) — > IR is differentiable with f convex, then for all 
m > 2 and all ti > (i = 1, 2, . . . , m), 



(m — l)(m — 2) 



l<«<j'<m 



<(m-2)j;/(t i )+/(x;t i ). (12) 



i=i \«=i 



Note that in (12), as well as in the sequel, the summation Y2i<i<j< m 
means that the sum is over all m(m — l)/2 pairs that satisfy 1 < i < 
j < n. We postpone the proof of the lemma. 

Proof of Theorem^ Suppose that a Steiner point of degree m+1 > 3 exists. 
We intend to show that m + 1 = 3. Note that we do not only consider the 
case m + 1 = 4, since it is a priori possible that degree 4 Steiner points don't 
exist, although degree 5 points exist (although we don't have any examples, 
and it seems highly unlikely). 

By Corollary [3] there exist vectors v\, . . . , v m and numbers £1, . . . ,t m > 
that satisfy d9]>, dTob and (fTTl). Square (TlOb: 



\ 2 



1=1 / 1=1 
in 



^|N|l + 2 ( v ^ v j) 
i=l i<*<i<™ 
m 

Y J (w{t l )f + 2 K«i> ( b y§»- ( 13 ) 

i=l i<*<i<™ 



Estimate (vi,Vj) by applying (11) to 7 = {i, j}: 

2 (vi,Vj) = \\vi + Vj\\l - \\vi\\l - \\vjWl 

< w(U + tj) 2 — w(ti) 2 — w(tj) 2 (again by 



13 



Sum this inequality over all pairs with 1 < i < j < m: 

2 («<» v i) ^ E M** + ^') 2 " ~ w{t i^ 

l<i<jr'<m l<i<j'<m 

in 

\<i<j<m i=l 

since each (u;(ij)) 2 is summed once for each of the m — 1 pairs in which i 
appears. Substitute this into (13): 

h(X>) ^ E Mii + i;)) 2 -(m-2)5>(^)) 2 - (14) 



l<i<j<m 

Apply Lemma [H] to f = w 2 : 
(m — l)(m — 2) 



i=l 



Mo)) 2 + 53 Wt i + ^-)) : 



l<i<J<m 



<{m-2)^(w(U)) 2 +[w(£ l U)) ■ ( 15 ) 



i=l 

(m— l)(rn— 2) / n ; n ^ 2 



1=1 



Combining (14) and (15), we obtain — 1, ' 2 v " t ^ (w(0)) < 0, which implies 
m + l<3. " □ 

Proof of Lemma^ We use induction on m > 2. The base case m = 2 is 
trivial. 

Assume now that m > 3 and that the lemma holds for m— 1, in particular 
that 



(m — 2)(m — 3) 



/(o)+ E + *i 



l<i<i<m— 1 



m—1 



m— 1 



<(m- 3) £/&) + /(!>)■ (16) 



i=l 



i=l 



Consider x := i m to be variable and i 1; . . . , i m -i fixed. Set T := 5Zi=i * 
and define 



S(a0 := (m - 2) E/fe) + /(!>)" E + 

j=l i=l l<i<j'<m 

m—1 

= (m - 2) ^ /(*<) + (m - 2)/(x) + /(T + z) 



i=i 



m— 1 



E E/(^ +a; )- 

l<i<7<m— 1 



i=l 
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We have to show that g{x) > 



(m-l)(m-2) 



/(0) for all x > 0. First of all, 



m—l 



g(0) = (m-2)J2 f{U) + {m- 2)/(0) + f(T) 



i=i 



m—l 



> 



l<i<j'<m. — 1 i=l 
m—l 

( m - 3) E /(*;) + /CO- E /(*i + *j) + ("»-2)/(o) 

i=l l<i<jr'<m— 1 

(m — 2)(m — 3) 



(m — l)(m — 2) 



./( ) + ( m _2)/(0) (by©) 
/(0). 



It is therefore sufficient to show that g'{x) > for all x > 0. We have 



m—l 



g'(x) = (m- 2) f{x) + f(T + x) - E /'(^ + 



i=i 



If T = then ij = for all i = 1, ... m — 1, and then is identically 0. 
We may therefore assume without loss of generality that T > 0. Write each 
x + L- as a convex combination of x and x + T: 



%+;e= ( 1-^ )x + ^(x + T). 



Since /' is convex, 



m—l 



m—l 



E m + x) < e ( 1 - 1 ) /'(^) + 1/'^ + r )) 



= (m-2)f(x) + f(x + T), 
which gives g'{x) > 0. This finishes the induction step and the proof. □ 



We now show that for each a G (1/2, 1) there exist Gilbert arborescences 
with Steiner points of degree 4 if the weight function is w(t) = d + ht a , with 
d, h > chosen appropriately. In the example all incoming tonnages will 
be equal. We first use Corollary [3] to formulate a result for general weight 
functions. 

Proposition 6. Let w(-) be a weight function that satisfies Conditions ([T|- 
Q. There exists an MGA with degree 4 in Euclidean 3-space with equal 
tonnages t\ = ti = t^ =: t and with weight function w(-) if, and only if 

3w(t) 2 + w(3t) 2 < 3w(2t) 2 . (17) 
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Proof. By Corollary [3j an MGA of degree 4 exists in Euclidean space with 
equal tonnages t\ = t 2 = t 3 =: t if, and only if, there exist three Euclidean 
vectors vx,v<2, v 3 such that 

IKH2 = \\v2W2 = IMI2 = w(t), (18) 
\\vi + V2W2, \\v2 + v 3 \\ 2 , \\vi + w 3 || 2 < w(2i), (19) 

\\v 1 +V2 + V 3 \\ 2 =w{3t). (20) 

(These vectors will of course span a space of dimension at most 3.) 



Square (20 ) and use ( 18 ) to obtain an expression for the sum of the three 



inner products: 

2(K v 2 ) + <u 2 , u 3 > + (vi,v 3 )) = w(3t) 2 - 3w(t) 2 . (21) 



Square (|19j) and use (18) to obtain an upper bound on each inner product 
{vi,Vj) (geometrically this is a lower bound on the angle between any two 
vectors): 

2 (vi,Vj) < w(2t) 2 - 2w(t) 2 , 



and substitute this into (21) to obtain (117 



Conversely, if we assume (17), we have to find three vectors that satisfy 
(18)-(20). Note that for each A £ [—1/2,1] there exist three unit vectors 
ui,u 2 , U3 £ M 3 such that the inner product of each pair equals A. The one 
extreme A = 1 corresponds to three equal vectors, and the other extreme 
A = —1/2 to three coplanar vectors such that any two are at an angle of 
120°. If we set 

w{3t) 2 1 



A 



6w(t) 2 2' 



then the vectors Vi := w(t)ui satisfy (|18[)-(20 ). It remains to show that this 
value of A really lies in the interval [—1/2, 1]. The lower bound A > —1/2 
holds trivially, while the upper bound A < 1 follows from < w(3t) < 3w(t), 
which in turn follows from non-negativity (Condition (fl|) and concavity 
(Condition Q) of the cost function. □ 

Corollary 7. For each a G (1/2,1) there exists d,h,t > and an MGA 
of degree 4 in Euclidean 3-space with cost function w(t) = d + ht a , and 
tonnages t\ = t 2 = t% := t. 

Proof. By choosing the unit of the weight function appropriately, we may 
assume without loss of generality that w(t) = D + t a . We may similarly 
assume that t = 1 , and then by Proposition [6] we only have to show that 
3w(l) 2 + w(3) 2 < 3w(2) 2 will hold for some value of D > 0, which is 

3(D + l) 2 + (D + 3 Q ) 2 < 3(D + 2 a ) 2 , 

or equivalently, 

D 2 + (6 + 2 • 3 a - 6 • 2 a )D + 3 + 3 2a - 3 ■ 2 2a < 0. 



16 



A sufficient condition for this to hold for some D > 0, is that the quadratic 
polynomial in D on the left has a positive root. For this to hold it is in turn 
sufficient that its constant coefficient is negative, i.e., that 

f(a) = 3 + 3 2a - 3 • 2 2a < for all a G (1/2, 1). 

However, it is easily checked that /(1/2) = /(l) = and that f"(a) > for 
all ol £ (1/2,1), so that / is convex on (1/2, 1). It follows that / is negative 
on (1/2, 1), which finishes the proof. □ 

6 Conclusion 

In this paper we have studied the problem of designing a minimum cost flow 
network interconnecting n sources and a single sink, each with known loca- 
tions and flows, in general finite-dimensional normed spaces. The network 
may contain other unprescribed nodes, known as Steiner points. For concave 
increasing cost functions, a minimum cost network of this sort has a tree 
topology, and hence can be called a Minimum Gilbert Arborescence (MGA). 
We have characterised the local topological structure of Steiner points in 
MGAs for linear weight functions, specifically showing that Steiner points 
necessarily have degree 3. 
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